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ABSTRACT 


An infinitely long, ring-stiffened, submerged, elastic cylinder 
having uniformly spaced elastic bulkheads is the structure considered. 
Loading is applied by a plane acoustic shock wave with front parallel 
to the cylinder axis. Dynamic pressure in the fluid is resolved into 
a free~field incident part and a scattered part. Structural response 
and scattered pressure in the surrounding fluid are found using finite 
element modeling of structure and fluid. Introduction of Fourier 
series makes the fluid region mathematically two-dimensional. A radia~ 
tion, or non~reflecting, condition at the outer boundary of the fluid 
ae is shown to give good results. A parametric study is made of 
effects of shock pulse rise time and duration on structural response. 
Results are presented as combinations of shock pressure and submergence 


pressure just sufficient to induce structural failure. 
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I. INTRODUCTION 


This investigation is the third in a sequence of finite element 
studies of fluid-structure interaction. The initial study was by 
Atchison [1] and the second by Newton and Atchison [2]. The system con- 
sidered, as in the earlier studies, is a submarine subjected to a shock 
wave resulting from an underwater explosion. The submarine is modeled 
as a ring-stiffened cylinder of infinite length with equally spaced 
elastic bulkheads. The shock wave is a plane acoustic wave with front 
parallel to the longitudinal axis of the cylinder. The structure-fluid 
equations are coupled by terms resulting from dynamic interaction at 


the interface boundary. 


A. PRIOR STUDIES 

It is shown by both authors in references [1] and [2] that the three- 
dimensional problem may be resolved into two-dimensional subproblems, 
acoustic wave propagation and transient response of a submerged struc- 
ture, using the finite element method. 

Acoustic wave propagation and spatial fluid discretization are pre- 
sented in [1]. The superposition theorem, wkich permits the rigid body 
pressure distribution to be obtained from a two-dimensional analysis is 
also presented in [1]. 

A refined structural model is used in [2]. In this paper, the cir- 
cumferential variations of fluid pressures and structure displacements 
are represented by Fourier expansions. Fluid pressures, structural de- 
flections, shell and ring stresses are tabulated for each harmonic 


through n=4. Extreme stresses resulting from superposition of harmonic 





contributions are given also. Two separate failure criteria are used 

to determine critical combinations of submergence pressure and shock 
pressure. The first of these utilizes limit analysis to evaluate combined 
effects of longitudinal and axial stress in the shell wall at a bulk- 
head. The second postulates that shell collapse will occur when the 


von Mises yield condition is satisfied away from the bulkhead. 


B. PRESENT STUDY 

In the present study the need for a separate solution of the acous- 
tic propagation is eliminated, thus effecting a significant reduction 
in computation. This improvement is made by dealing separately with 
the incident shock pressure - an entity which is unaffected by the pres- 
ence of the structure. The remainder of the dynamic fluid pressure 
(called the "scattered" pressure) is represented in the finite element 
modeling of fluid-structure response. As in [2], the azimuthal varia-~ 
tion of fluid pressure and structural displacement are represented by 
trigonometric series so that the fluid region modeled is mathematically 
two-dimensional. 

A further improvement introduced here is the imposition of a radia-~ 
tion boundary condition at the cylindric outer boundary of the fluid. 
Although this condition is necessarily approximate, it is found to re- 
duce significantly the effects of boundary reflections on structural 
response. 

A systematic study is made of the effects of shock pulse duration 
and rise time on structural response. Results are presented in the form 
of critical combinations of shock pressure and submergence pressure se 


sufficient to induce structural failure. 





II, WAVE PROPAGATION 


A. PRESSURE RESOLUTION 
A submerged structure subjected to an underwater shock wave experi- 
ences,at a given time, a total dynamic pressure which can be considered 


to have three components 


where 
p is the total dynamic pressure in the fluid, 
p. is the pressure due to the incident wave, considering 
the structure absent, 
Pp is the pressure due to the reflected wave, consider- 
ing no structural deflection, 
is the pressure radiated by the structure due to 


Da 


structural deflection. 


There are two useful combinations of these pressures for purposes 


of analysis. Consider first 


Prig 5 Pi , Pr (2) 


where ae represents the total dynamic pressure field with a rigid 


structure. Next consider 


fe wep Py (3) 


where Ps represents the scattered pressure, which includes both re- 


flected and radiated parts. 


LO 





Instead of considering the total dynamic pressure to be the sum of 


Di... and p-_ as in [1] and [2], the resolution 
rig Yr 


p= par Pp (4) 


will be used. This resolution has been employed widely in similar prob- 
lems, e.g., by Mindlin and Bleich [3], Geers [4] and Everstine [5]. 
The discretized equations for the structure and fluid given in [1] 


and [2] are 


Mé+K6=5Lp (5) 
Qp+Dpt+Hp=-pL 6 (6) 


where § is a column vector of structural nodal displacements and p 

is a column vector of fluid nodal pressures. M, K, Q, D and H are sym- 

Metric coefficient matrices, L is a matrix describing coupling at the 

structure-fluid interface, and p is the fluid density. Superior dots 

denote time derivatives and superscript T denotes transposition. 
Considering first the incident pressure only, a time-dependent column 

vector f£ is added to correct for the unsatisfied boundary condition 


at the fluid-structure interface. It is required that 


[He 


Op, +Dp, +Hp, =f (7) 


1 1 


Considering next the scattered pressure, it is required that 


Ito : 


° 2 T a 
QP. tDepet+ Hype pb 6 - = (8) 


S == 
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The equation for the structure (5) becomes 


{Or : 


mMO+K6=Lp, +p. (9) 


Therefore, equations (7), (8) and (9) are the new equations for the com- 
plete elasto-hydrodynamic problem formulated by Atchison [1] and Newton 
and Atchison [2]. 

Rather than use equation (7) to find f£, the relation between the in- 
cident pressure and the fluid particle velocity normal to the boundary 
may be used. 

Consider a plane wave traveling in the negative x-direction. The 


free-field incident pressure is given by 


Pp, = g(x + ct) {10) 


where c is the acoustic velocity of propagation. Fig. 1 represents. 


the wave front. 


x+ct 


Fig. 1. Wave front. 


Applying Newton's second law to a fluid particle, it can be shown 


that the fluid acceleration is given by 


© coco tee (11) 





However 


= at = 5 
Pi oe. 2 oC 
Therefore, in (11) 
* 1 e 
Vv a pc se (12) 


Formula (12) gives the relation between the incident pressure and the 
fluid particle velocity. The complete problem now is reduced to solv- 
ing the set of linear equations given by (8) and (9) using the result 
established by equation (12). 

Fourier representations are employed in deriving the discretized 
counterparts of equations (8) and (9) for each pressure harmonic and 
displacement harmonic. Because both the structure and the fluid region 
are axisymmetric, there is no interharmonic coupling. 


Therefore, for the structure displacement, let 
6 = 56 cos n@ (13) 
n 
and, for the fluid pressure, let 
p reo cos n@ (14) 


where n=0,1,2,... and summation with respect ton is implied. 

For each Fourier order n there is a separate pair of equations (8) 
and (9) governing the Soe ee ead pair Oe and ee of displacement 
and pressure vectors. The coefficient matrices of these equations are 


order-dependent. Formulas for these matrices are given in Appendix A. 


B. RADIATION BOUNDARY CONDITION 
Finite element modeling of infinite regions always requires some 


Sstrategem to avoid significant errors due to substitution of a finite 
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region. In most static problems and some steady-state dynamic problems 
it suffices to place the boundaries "far enough" from the region of 
interest. This scheme is also applicable to the present problem featur- 
ing wave propagation in the fluid. Here the requirement of "far enough" 
is satisfied if waves traveling from structure to outer boundary and 
back to the structure arrive after the occurrence of maximum stresses in 
the structure. Meeting this requirement led, in the present applica- 
tion, to a finite element model in which over 80 percent of the degrees 
of freedom represent fluid pressures. 

An alternate means for simulating an infinite region utilizes a 
radiation, or non-reflecting, condition at the outer boundaries. Such 
a condition was first proposed by Zienkiewicz and Newton [6] in a form 
which is exact for plane waves normally incident upon the boundary. 
Chenault successfully employed the condition in a study of two-dimensional 
added mass and damping [7]. Dean constructed and successfully demonstra- 
ted an extension to spherically spreading waves [8]. For both plane 
and spherical waves the corresponding non-reflecting conditions are 
exact, but no such result is possible for two-dimensional (cylindric) 
spreading. Despite the necessity for approximation, Bai constructed and 
used a radiation condition for two dimensions [9]. An extension of Bai's 


scheme is developed in the section entitled "Two-Dimensional Analysis". 


C. ONE-DIMENSIONAL ANALYSIS 

A one-dimensional analysis using the new pressure resolution has been 
done in order to compare new results with those previously obtained by 
Atchison [1]. 

The region considered is a semi-infinite fluid strip of unit cross 


sectional area. Through this region a shock wave of known shape is 


14 





propagating. The semi-infinite strip is modeled with a one degree of 
freedom structure at one end and the radiation boundary condition is 
applied at the other end. 

The results obtained were compared with those obtained by using 
CSMP (continuous system modeling program). Good agreement was found. 
Accordingly, the Houbolt,integration method [10], the linear element, and 
the restrictions for node spacing and time step as used in [1] and [2] 
were found to give excellent results when used with the new pressure 


resolution. 


D. TWO-DIMENSIONAL ANALYSIS 
1. Boundary Conditions 
A non-reflecting or radiation boundary condition is presented 
and used with this fluid-structure investigation. 
If it is wished to let a plane wave traveling in ete positive x- 
direction pass out of the region being investigated, the condition re- 


quired is 


dp 


LE p (15) 


Qtr 


In dealing with two-dimensional problems, the axisymmetric ap- 
proximation proposed by Bai [9] is adapted. The (scattered) pressure 
distribution for a wave moving in the positive r-direction is assumed to 


be described by 


wel /2 


p 2 g, (e-ct) cos n@ (16) 


Ss 


where, as before, there is summation with respect to index n. 


JUS: 





The boundary condition imposes a requirement on the normal deri- 
vative of the pressure at the boundary. The partial derivative is cal- 


culated from equation (16) as 


dP 
S| db -3/2 = -1/2 " 
a oe g_(x-ct) cos nO +r gi (x ct) cos n@ (17) 
Also 
op. 
2s ~ =lf2 » a 
a gi (x ct) cos n@(-c) 
or : 
meee L/2 : 
a xr gi (x ct) cos né@ 


Therefore, equation (17) may be rewritten as 


dp p p 
a agin (18) 
or 2x Cc 


This formula will be introduced as the radiation boundary condition for 
the solution of the interaction problem. Data concerning the usefulness 
of this approximation are given in the section titled Results. 
2. Fourier Analysis 

It is required, in order to solve equations (8) and (9), to know 
both the incident pressure and the normal component of the fluid acce- 
leration at the structure-fluid interface. 

The pressure-time history of the cubic ramp-box wave used may be 
represented as shown in Fig. 2, where TY is the pulse rise time, T 


d 
is the pulse duration and To 1s the subsidence time. 
For a given incident pressure wave, characterized by rise time, 


pulse duration, subsidence time, acoustic velocity and fluid density, it 


is required to find Fourier components for the incident pressure and 
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Te 79 Ts : 


Fig. 2. Pressure-time history of a cubic ramp~-box wave. 


fluid acceleration at the shell surface using a chosen time step h. 
These Fourier components are needed for as many time steps as are used 
in solving the fluid-structure equations. 


Basic parameters for the Fourier analysis are shown in Fig. 3. 





Fig. 3. Parameters used in Fourier analysis. 


Ly 





For the cubic wave front, pressure and pressure derivative formulas are 


derived by Atchison in [1] 


2 2. * 
36 =- 26°7p ’ = < in 

Pp; = (19) 

. < < 

p ? tS eS aa) 
6p* 2 
ic -t*) : ec 
T r 

o i a 

P, = (20) 
0 , Cee eee AT et) 


where € = = and p* is the shock pressure. 

The eetagnce is described by the same kind of cubic transition 
as the rise. 

Relations for the normal (perpendicular) component ve of the 
fluid acceleration can be derived uSing equations (12) and (19). Thus 


6 p* “2 
pe T, (Z GC) cosé , co TL 


p> (21) 
0 


< 
pees baer TS) 


Therefore, equations (19) and (21) are used as basic formulation 
for the Fourier analysis. 

Fourier components for the incident pressure and fluid accelera- 
tion are calculated for each time step using 24 intervals over the 
cylinder from 9 = 0° to @ = 180°. It is not necessary to include values 


of 8 from 180° to 360° because of symmetry about 6 = 0°. 
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For the incident pressure coefficients 


N 
re) N J 15 
J=0 
n= 41,2,3,4 
N 
eS 2 w_p. cos n§@ 
n N Tee J 
J 
J=0 
where We = Wis = 1/2 and all other We = 1. 


For the fluid normal acceleration components 


N 
ee — N = 24 (23) 
O N aj Ps 

> eee Z 

N 
B “— : fa) 
n N Ws i cos n J 

J=0 


Required values of P, and x are. found from equations (19) 


and (21). 


New Fourier resolutions are required when any characteristic 


time of the shock pulse is changed. 


ane 





III. STRUCTURAL RESPONSE 


A. STRUCTURAL MODEL 

A refined model is used for the structure. A total of eight struc- 
tural elements and 27 structural degrees of freedom per Fourier component 
represent the final mesh and it is identical with the one utilized by 
Newton and Atchison [2]. 

As described in [2], the structure is a ring-stiffened elastic shell 
of infinite length with elastic bulkheads at spacing 2L. Fig. 4 shows 
a cross-section of the repeating length of structure and the associated 


fluid region. 


N 
- 
o* 


t 


Soe ee SS OS Se 


Fig. 4. Region represented by finite element model. 


It 1s assumed that the plane of a bulkhead, z=0, and the plane mid- 
way between bulkheads, z=L, are planes of symmetry for the structure 
and the response displacements. 

As in reference [2], the effects of the stiffening rings are in- 


Cluded by treating the shell as orthotropic. Stress calculations take 
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into account the geometry of the rings. The effective shell thickness 
in the circumferential direction includes a contribution by the rings. 

The bulkheads furnish elastic restraint against radial displacement 
of the shell. This is accomplished by treating the bulkhead as a 
diaphragm of uniform thickness. 

As developed in reference [2] a geometric stiffness matrix modifies 
the structural stiffness. It considers the effects of the membrane 
compressive stresses of the axisymmetric mode acting to amplify deflec- 
tions in the higher harmonics (n>1). 

Further information concerning the structural model is given in 
reference [2]. Details concerning application of Fourier series to 


analysis of axisymmetric shells are given by Grafton and Strome [11]. 


B. FLUID MODEL 

As indicated in Fig. 4 the length of the fluid veeion is one-half 
of the bulkhead spacing, its inner radius is the shell radius "a", and 
the outer radius is Xr. 

Rectangular elements having corner nodes and linear shape function 
are used. Along the fluid-structure interface a fluid node is located 
at each structural node. 

The radiation boundary condition as developed previously is applied 
at the outer radius roe For the purpose of studying the effectiveness 
of the radiation boundary condition, the same combination of pulse para- 
meters and submergence pressure has been employed using a=5 meters and 
a sequence of values of ry from 8.15 meters to 17.6 meters. 

In finding critical combinations of shock and submergence pressure 


to induce structural failure two different meshes are utilized. In each 
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case the fluid region has inner radius a = 5 meters and outer radius — 
= 16.025 meters, The first mesh has 90 fluid elements and 120 fluid 
nodes and each element has a thickness (radial) of 0.525 meters. The 
second has 174 fluid elements and 225 fluid nodes with element thickness 
0.2625 meters. The latter mesh is used when the rise time of the wave 
front is small in order to meet the restriction given by Atchison [1] as 


Gu 
r 


s < oe (24) 


where s is the node spacing. Two different time steps are employed 
depending on the size of the mesh. For the finer mesh a time step of 
O.1 milliseconds is used. For the coarser mesh the time step is 0.3 


milliseconds. 


C. SOLUTION TECHNIQUE 

The solution process iS very similar to the one used by Newton and 
Atchison [2]. 

There iS a separate solution for each of the five harmonics considered 
(from n=0 to n=4). Once the solution for the first harmonic is completed, 
the solution for the second starts and the other harmonics follow in 
order. 

Fourier components for the incident pressure and normal fluid 
acceleration are introduced as input for each solution. 

For any harmonic the procedure employed first advances 6 by one time 
Step using equation (9). The value of p, at the forward step is obtained 
by parabolic extrapolation. With 6 and 8 at the forward point, equa- 
tion (8) is then used to make the forward step for Poe Successive 
iterations are performed using the latest values of p. and §-on the right 


hand sides of equations (9) and (8). 
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Detailed information about this solution technigue is given in refer- 


ence [2]. 


D. 


MODEL PARAMETERS 


The following parameters are used for structural response calcula- 


tions. 


Shell radius a 5 meters 


Distance between bulkheads : 16 meters 


Shell thickness 0.05 meters 


Circumferential average thickness : 0.07 meters 
Bulkhead thickness : 0.02 meters 


Structural material Steel 


Yield stress for steel 700 Megapascals 


Shell flexural rigidity ratio Dp/D_, : 100 


Shell twisting rididity ratio Do, /D_: Oso 


FAILURE CRITERIA 


Failure criteria used by Newton and Atchison are employed. The 


following excerpt from reference [2] defines the basis. 


The largest calculated elastic stresses are found to be 
in the shell at the bulkhead. These stresses result from 
Superposition of longitudinal bending and axial membrane 
compression. Since yielding at this location would not pre- 
cipitate shell collapse, combinations of bending moment and 
axial force corresponding to failure are based on limit 
analysis (full yielding through thickness - part compression, 
remainder tension). 

A different situation exists at appreciable distances 
from the bulkhead. Here it is postulated that yielding will 
induce shell collapse. At such locations an effective uni- 
axial stress 0 is calculated using the von Mises yield condi- 
tion. Failure is deemed to occur when the utilization ratio 
R= eer reaches unity. 
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IV. RESULTS 


A. INTRODUCTION 

Results reported here concern two separate investigations. The first 
of these is a limited examination of the effectiveness of the radiation 
boundary condition. The second, and principal, study is a systematic 
examination of the effects of shock pulse rise time and duration on 


structural failure. 


B. EFFECTIVENESS OF RADIATION BOUNDARY 

The effectiveness of the radiation boundary condition is examined 
by repeated solutions of the same problem using successive reductions in 
the outside radius ro of the fluid region. Table I gives the correspond- 
ing utilization ratios for both failure modes. The first line of the — 
table provides the comparison standard, since the extreme values of the 
utilization ratios occur at 15.6 milliseconds and the time for an 
acoustic wave to travel from shell to the outer boundary and back is 18 
milliseconds. The second line (ro = 16.025 meters) also represents a 
case where the outer boundary is "far enough" from the shell so that 
the utilization :atios are unaltered. Succeeding radius reductions 
do affect the utilization ratios, but even for the smallest region (xO 
= 8.15 meters) the ratios change by less than 3 percent. The accompany- 
ing reduction in the number of fluid nodes is dramatic. In reducing 
from ry = 16.025 meters to ro = 8.15 meters, the number of fluid nodes 


goes from 120 to 45. 
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TABLE I. Effect of Radiation Boundary Placement on Utilization Ratios 


OUTSIDE RADIUS TRAVEL TIME RELATIVE UTILIZATION RATIOS 
(METERS ) (MS) LIMIT ANALYSIS VON MISES 





76 1.000 
16.025 1.000 
14.45 0.998 
12.875 0.997 
Jl gas 0.976 
I. 725 0.974 
8.15 0.971 
Notes: 


1. Shell radius: a= 5 meters. 
2. Cubic ramp-step wave with ctT,/a = 0.97. 
3. Submergence pressure = 4.0 Megapascals, 
shock pressure = 4.5 Megapascals. 
4. All maxima occur between 14.3 and 15.7 milliseconds after onset. 
5. Travel time is time for acoustic wave to travel from shell to 
outer radius and back to shell (c = 1.4 meter/millisecond). 


TABLE II. Effect of Fixed Outer Boundary Placement on 
Utilization Ratios 


OUTSIDE RADIUS TRAVEL TIME RELATIVE UTILIZATION RATIOS 
(METERS ) (MS) LIMIT ANALYSIS VON MISES 


17.6 LB ..0 1.000 1.000 
9.725 6.8 On993 0.949 
8.15 4.5 Of 792 0.947 


Notes: Notes for Table I apply here. 


M8. 





Table II exhibits the effect of reducing the outer radius of the 
region without imposing the radiation condition at that boundary. It 
is apparent that the resulting degradation of accuracy is more severe 


than with the radiation condition. 


C. INFLUENCE OF PULSE RISE TIME AND DURATION ON STRUCTURE FAILURE 
Considered here is the principal application of the improved analy- 
sis techniques introduced in this thesis. In reference [2] results 
were obtained for structural response to a cubic ramp-step shock wave. 
A single rise time ct /a = 0.97 was considered and the duration was 
effectively infinite. A large family of cubic ramp~box pulses is 
treated here. The rise time TY and effective duration T = > T +t, + 
i, are systematically varied. (Note that the impulse of the shock 


PAS 
pulse is p*T. See Fig. 5.) 








Te Ta Ts 


Fig. 5. Effective pulse duration T. 


The pulse durations studied here are cT/a = 2,4, and 8. For each 
duration, four separate rise times: ct /a = 0.1,0.2,0.4, and 0.8 are 
considered. In every instance the fluid region extends from a = 5 


meters to Ly = 16.025 meters. For all combinations with ct /a > O.l a 
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fluid mesh with 90 fluid elements and 120 fluid nodes is used "coarse" 
mesh). The thickness (radial) of each element is 0.525 meters. With 
this mesh a time step h = 0.3 millisecond is used and the subsidence 
time is cT./a = 0.2. For the rise time ct /a = 0.1 a "fine" mesh is 


introduced with 174 fluid elements, 225 fluid nodes and a 0.2625 meter 


element thickness. With this mesh a time step h 0.1 millisecond is 


used and the subsidence time is reduced to cT./a Dede 

For both meshes there are eight shell elements with lengths varying 
from 0.5 meters near the bulkhead to 2.0 meters midway between bulkheads. 
The number of structural degrees of freedom is 27 for each harmonic. 
Along the fluid~structure interface there is a fluid node coinciding 
with each structural node. 

Results are presented graphically in the form of failure boundaries 
on the submergence pressure versus shock pressure plane in Figs. 6 
through 17. Numerical results used to define these See ies are tabu- 
lated in Appendix B. 

Fig. 7 presents two sets of failure boundaries for ct /a = 0.2, 
cT/fa = 2. For each boundary the narrow line is found using the fine 
mesh and the wide line is obtained from the coarse mesh. 

When the effective pulse duration is small (cT/a = 2), changes in 
pulse rise time affect considerably the failure boundary. Increase in 
effective pulse duration (cT/a > 4) produces failure boundaries that are 


not sensitive to changes in rise time. In fact, for cT/a = 8 the same 


O:d. Tt 1s 


values are obtained for all rise times greater than ct /a 
considered that the small difference observed for ct, /a = 0.1 is largely 
due to the mesh change. For durations greater than cT/a = 4 results 


approach those of Fig. 16, regardless of rise time. These results are 
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based on ct /a = 0.97 and cT/fa = ~, the pulse parameters used by Newton 
and Atchison [2]. The boundaries found here agree closely with those 
of reference [2], but it is noted that labels on ordinate and abscissa 


have been mistakenly interchanged on Fig. 10 of [2]. 
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V. CONCLUSIONS 


It is concluded that the separation of the total dynamic pressure 
acting on a structure into incident and scattered pressures has im- 
proved the analysis of the fluid-structure interaction problem investi- 
gated in this thesis by avoiding the propagation solution. The results 
obtained also confirm that the imposition of the radiation boundary 
condition atthe cylindric outer boundary affords a further improvement. 
For a representative structure, the effects on failure of pulse rise 


time and pulse duration have been established. 
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APPENDIX A. FLUID COEFFICIENT MATRICES 


1. Separation of fluid equations 


The governing equation for the scattered pressure is 
Qp+DB+Hp=PL S8-£ (8) 


where p is the vector of nodal scattered pressures. It 1s assumed 
that the (scattered) pressure p ata point (r, 9, z) of the fluid is 


given by 
p= N. P,,, cos nd (Al) 


where N. , a function of r and z=, is the shape function associated 
with fluid node i and Pin is the nth Fourier component of pressure at 
that node. Summation on indices i andn is implied. 

This leads to a separate equation like equation 8 for each Fourier com- 


ponent of vector p. The typical equation may be written as 


2,8, +2, 8, +H 2,8 8S (aa 


where the summation convention does not apply. Equation A2 assumes that 


Fourier resolution has also been applied to the structural nodal displace- 


Ments and that the normal displacement w of the shell is given by 


w= NN," 6... cos no (A3) 
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where N,! , a function of z=, is the shape function for structural 
node i and On is the nth Fourier component of displacement at this 


node. As in equation (Al) the summation on i andn applies. 


2. Formulas for matrix elements 

Coefficient matrices are found by calculating contributions at the 
element level and assembling by addition. Presented in this section 
are the element-level contributions based on the convention that a lower 
case letter represents an element of the matrix denoted by the corre- 
Sponding capital letter. Thus, aan is the element in row i, colum jj, 
of matrix Q- 

In each of the following formulas a symbol e , the Neumann factor, 
appears. It is defined to have the value 1 for n=0 and the value 2 for 
nN any positive integer. In the two following formulas the integrations 


with respect to r and z extend over the range for the individual fluid 








element. 
20 
a. oS J{fN. N. xv ar dz (Ad) 
ijn : ea 
EC 
n 
OT ON. ON. ON. ON. a2 
h.. = So tye et ett ee.) rar az = (a5) 
ijn om dr or dz 02 2 ass 


No summation on n is intended in equation (A5). 

The following two formulas concern the boundary condition at the 
structure-fluid interface (r=a). The only fluid nodes having non-zero 
contributions are those in contact with the shell. Integrations with 
respect to z take place along a structural element and the matching 


fluid element. 
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o. = fo on.tn. dz (A6) 
ijn ¢€ » 4 

n 
eee! ae ee (A7) 
jn om j 


The component Pan , when multiplied by BL as defined in equation (23), 
gives the jth component of vector - fn : 

The final two formulas represent contributions due to the radiation 
boundary condition applied at the outer boundary (r = ra): Integration 
with respect to Zz is carried out along that boundary. Only fluid nodes 


on that boundary are included 





amr 4 
cl f WN. N. dz (A8) 
ijn cee All g 
T 
Ah.. =——~— Jf N.N, dz (A9) 
aya 0) 


The contributions Ah.. are assembled into H . 
ijn —n 
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